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Abstract 

The goal of this paper is two-fold. First, based on the interpretation of a quantum tight-binding 
model in terms of a classical Hamiltonian map, we consider the Anderson localization (AL) problem 
as the Fermi-Pasta-Ulam (FPU) effect in a modified dynamical system containing both stable and 
unstable (inverted) modes. Delocalized states in the AL are analogous to the stable quasi-periodic 
motion in FPU; whereas localized states are analogous to thermalization, respectively. The second 
aim is to use the classical Hamilton map for a simplified derivation of exact equations for the 
localization operator H{z). The letter was presented earlier [J.Phys.: Condens. Matter 14 (2002) 
13777] treating the AL as a generalized diffusion in a dynamical system. We demonstrate that 
counter-intuitive results of our studies of the AL are similar to the FPU counter-intuitivity. 
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I. INTRODUCTION 



Half a century ago two celebrated papers were published temporally close to each other 
which gave birth to two fundamental directions of theoretical physics. In 1955 the Fermi- 
Pasta-Ulam (hereafter, FPU) paradox was formulated 1, ^ which suggested the nonequipar- 
tition of energy among normal modes of an anharmonic atomic chain. This phenomenon 
is closely connected with the problems of ergocidity, integrability, chaos and stability of 
motion [2, jj. A few years later, in 1958, Anderson [J] suggested the possibility of electron 
localization (AL) in a random system, provided that the disorder is sufficiently large. This 
idea is one of the foundations for the understanding the electronic properties of disordered 
systems jg]. From a more general point of view, this idea implies the absence of wave dif- 
fusion in a random medium as a universal feature of stochastic processes. Unlike the AL 
problem, where the stochasticity was explicitly introduced into consideration through ran- 
dom potentials, the FPU considers the stochasticity as a non-trivial effect in the dynamics 
of nonlinear systems, since a strong chaos behaviour can be observed even in a system with 
several degrees of freedom. 

Despite the fact that the FPU paradox can be reformulated analogously to the AL prob- 
lem in terms of waves (ineraction of normal modes which are characterized by their wave 
number k), their conceptual similarity has not been noticed so far. In fact, localization, even 
when mentioned in the FPU literature, is not associated with disorder; usually this means 
that the energy initially placed in a low-frequency normal mode k$ of the linear problem 
stay almost completely locked within a few modes neighboring the ko mode. The local- 
ization in /c-space of normal modes [6| or analogous energy localization in the FPU chain 
j?! are typically considered here. Sometimes the AL and FPU are mentioned together, but 
as independent and even competing processes. In particular, energy transport in binary 
isotropically disorded nonlinear FPU chains was considered [8| with the competition be- 
tween localization (a disorder effect) and mode transitions (a nonlinearity effect). Notice 
also that the close connection between the nonlinear dynamics and the AL was established 
for nonlinear systems with a much smaller number of freedom degrees than in the FPU 
problem. Thus, it was shown that the quantum kicked rotor model (9), [lo| can be mapped 
onto the AL model. On the other hand, the quantum kicked rotor serves as the starting 
point of the systematic analysis of the quantum dynamics of classically chaotic dynamical 
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systems. 

Our purpose is to establish a close connection between the two fundamental problems. 
That is, the AL is nothing else but the FPU effect in a modified dynamical system with 
interacting normal modes. Besides, the delocalized states in the AL problem are analogous 
to the stable quasi-periodic motion (recurrence in the FPU problem) and, respectively, the 
localized states are similar to the thermalization (motion instability) in the FPU. Our FPU 
modification includes: (i) non-trivial change of the mode ensemble. In the FPU without 
interactions all normal modes are stable. In contrast, in the AL both stable and unstable 
(inverted) modes also exist, (ii) a modified mode interaction: the stochasticity in the AL 
is introduced directly, through random forces linear in coordinates and rather indirectly, 
through nonlinear terms in the dynamical equations. 

A critical comparison of the AL problem and the FPU paradox is the more useful since the 
former turned out to be also a paradox. Indeed, in recent years, the conflicting situation was 
established here: (a) experimental results contradict the generally-accepted theory 

HQ, 

whereas (b) the analytical theory contradicts numerical simulations 0. 



131 ] . Despite such a 



clear conflict, the results of our exact analytica 



unexpected in the Anderson community 



17. 
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theory [jj], Q, Q are considered as highly 
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20] since they contradict the generally- 



accepted results of the scaling theory [21] and numerical modeling [121 ] ■ In particular, doubts 
are expressed about our conclusion on the existence of the metal-insulator transition in 
the two dimensional (2D) disordered system of noninteracting electrons (which does not 
contradict real experimental data 0, [llj]) and that the Anderson transition is of the first- 
order (localized and conducting state co-exist). The latter statement is a clear example of 
a counter-intuitive prediction. Appearance of highly unexpected results is typical for such 
counter-intuitive (paradox) problems such as the AL and FPU. 

In this paper, we employ the main results of our analytical method Q, Q, Q and treat 
the AL as generalized diffusion in a dynamical system. The random forces impose random 
walk amplitudes (the dynamics is bounded in the phase space for unperturbed system) which 
can lead to the diffusion divergence. The diffusion problem with random forces linear in the 
coordinates can be exactly solved. Moreover, the diffusion concept permits to connect the AL 
and FPU problems. As is well-known [lOj, the equations for nonlinear dynamical systems 
can under certain conditions describe pseudo-random walks which leads to the diffusion 
behaviour and diffusion increase in mean energy (i.e. divergence). 
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We shall employ also the interpretation of the Schrodin ger equation (quantum tight 



binding model) in terms of the classical Hamiltonian map 22|, [23|, |24J . As a result, the prob- 
lem can be reformulated in terms of interacting mode dynamics which opens the opportu nity 
for the detailed comparison of AL with FPU. It should be stressed that ideas 



allow to simplify considerably the mathematical formalism used earlier 



flfl 




1611 . This is 



important in the light of the recent criticism [l?], that the engineering language (signal 
theory) of the mentioned formalism (input and output signals, filter function, etc.) is new 
for the AL community. In this paper we suggest a compact derivation of the main feature 
of the disordered system - the localization operator H(z) 16]. 

The structure of the paper is as follows. In Section III Al we explain how the one- 
dimensional tight-binding model with diagonal disorder can be presented in terms of the 
classical two-dimensional Hamiltonian map for normal or inverted oscillators. In Section 
III Bl we present the equations for arbitrary dimensions. We show that delocalized states, 
in general, correspond to the statistically bound trajectories, whereas localized states to 
unbound trajectories, respectively. The trajectory type depends on the excited mode (nor- 
mal or inverted). We show in Section [II CI that statistically unbound trajectories can be 
treated in terms of a generalized diffusion, with exactly predicable properties. As a result, 
we arrive in Section III Dl at the definition of the localization operator H(z). In Section [TTT1 
a possible comparison of the AL and FPU problems is discussed. It is shown that the AL 



treatment in terms of the classical Hamiltonian map resu 



ts in a paradox demonstrating the 



counter-intuitive nature of this problem. The results [14 
for the FPU and their detailed similarities are analyzed 



161 ] are compared with those 



II. ANDERSON LOCALIZATION AND CLASSICAL HAMILTONIAN MAP 



A. Cauchy problem and classical Hamiltonian map for the one-dimensional case 



As is well known, in order to determine the Lyapunov exponent 7 (which is the inverse 
of the localization length, £ = I/7) and the phase diagram (the areas of the localized 
and delocalized states), the Cauchy problem with fixed initial conditions has to be solved 



251 ] . For illustration, the ID Schrodinger equation 



(1) 
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with random potentials e n can be presented as a recursive relation 

V> n +1 = Elp n - - ^nlpn- (2) 

The treatment of one of the spatial coordinates as a temporal variable (discrete time n) is a 



standard approach in chaos theory 26J which opens the way to the dynamical interpretation. 



Taking into account the so-called causality principle [14], ll6j| , the latter equation permits an 
exact stochastic analysis. Indeed, it is easy to see that in the Cauchy problem (with fixed 
initial conditions for ipi and ipo), eq.(T5]), ip 2 is a function of ei, ^3 is a function of £2,^1, etc 
(a causality). That is both amplitudes ip n and ip n -\ on the rhs of eq.([2]) are statistically 
independent of e n and can be averaged separately (causality principle): 

^ n+1 ) = (E-(e n ))(ij n )-(4> n „ 1 ), (3) 

(r n+l } = (E 2 -2E(e n } + (el))m- (4) 
-2(£-(£„,))(VvAn-i> + (^-i>- 

The causality principle can be used only (a) for the recursive relation for the Cauchy problem 
and (b) when on-site potentials are independently distributed ((e n e n >) = a 2 5 n ,ri) but not for 
the Dirichlet problem, correlated potentials, etc. 

To study the origin of localized/delocalized states, we use a simple approach 
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based on the interpretation of a quantum tight-binding model, eq.([2]), with diagonal disorder 
in terms of the classical two-dimensional Hamiltonian map. The difference equation (T5]) is 
reduced to a discrete transform with a simple physical interpretation. Let us write the 
second-order equation as a set of two first-order equations. Assume q n = ip n and p n = 
■0n+i — Vv for E > 0; whereas q n = —ip n and p n = —{jp n+ i — i[) n ) as E < 0. The obtained 
equation set reads 

Pn+i =Pn~ u 2 q n - e n q n , (5) 
q n +i = q„.+Pn+i, (6) 

where uj 2 = 2 — \E\ . 

As is shown [26j, the discrete transform (jSJ), can be connected with the equivalent 
differential equation of the Hamilton dynamics with the Hamiltonian (kicked oscillator) 

H=£ + ^ + ^ £'<«-»*>■ (7) 



It defines the system with the unperturbed Hamiltonian for an oscillator affected by a 
periodic sequence of a kicks (5-pulses) with the period At. Discrete transform arises when 
magnitudes of coordinate and impulse are considered with a discrete time increment At and 
correspondingly with a discrete time index n (discrete time). 

Dependent on the uo 2 sign, there exist two cases. For \E\ < 2, corresponding to the 
solutions inside the band in the unperturbed system, we obtain a normal (stable) oscillator 
with oj 2 > 0. On the contrary, the energies outside the band, \E\ > 2, are associated 
with inverted (unstable) oscillators with uj 2 < 0. Without disorder e n = and delocalized 
quantum states correspond to a normal oscillator and trajectories bound in the classical (p, q) 



space as n — > oo 



231 ] . Simple characteristics of the dynamical system q 2 or p 2 n are bound, 



respectively. Contrary, an inverted oscillator describes non-physical solutions outside the 
band which are now unbound in the classical phase space; the q 2 magnitude is divergent as 
n — > oo. 

When disorder is introduced, the situation changes qualitatively. Random kicks for a 
normal oscillator lead to the random amplitude walks. The oscillatory motion remains since 
the average (q n ) remains bound as n — > oo. For a random amplitude walk long trajectories 
in the classical phase space (p,q) are possible which can be treated as diffusion 16|. This 
is characterized by a typical parameter divergence: (q 2 ) — > oo as n — > oo. Since these 
trajectories correspond in ID to the localized states, these can be considered as statistically 
unbound. In the case of the inverted oscillator the amplitude increases exponentially (in the 
model with Hamiltonian eq.Q) between successive kicks, but the force linear in coordinates 
is able to change a coordinate sign. 

Use of the diffusion terminology is quite justified here. Indeed, to detect the diffusion, 
it is sufficient to demonstrate the divergence of the second moment of the amplitude q n 
and to establish its time- dependence, the function f(n) in (q 2 ) = f(n). The divergence 
of the second moment defines the conditions of the diffusion appearance. For a normal 
diffusion the mean square displacement is linear in time, (q 2 ) oc n. The notion of an 
anomalous diffusion 27|, |28[ derives from the fact that the mean square displacement may 
be anomalously diffusive, (q 2 ) oc n a (a ^ 1), i.e. nonlinear in time (power-law divergence). 
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25|: (q 2 n ) oc exp(2 7 n) 



The quantity (q 2 ) = (ip 2 ) in eq.(j2D was calculated analytically 
with 7 > for an arbitrary E value, provided a > 0. In this case the generalized diffusion 
(exponential divergence) takes place. The appearance of the localization in the approach 
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based on eq. <$fy is equivalent to the appearance of diffusion. Respectively, the well-known 
statement that in one dimension all states are localized at any level of disorder, is equivalent 
to the statement on the diffusion character of all solutions of eq.© in ID for a > 0. 



B. Classical Hamiltonian map for D-dimensional case 

The above-discussed statement of the problem can be naturally generalized for an arbi 



trary dimension D; this idea was mentioned but not realized in Ref. |24|. The phase diagram 
of the system with metal-insulator transition should be obtained in the thermodynamical 
limit (the infinite system). Let us consider the semi-infinite system, or an infinite system 
with a boundary, where the index n = m£> > 0, but all rrij G (—00,00), j = l,2,...,p, 
with p = D — 1. We combine indices in the form of a vector m = {m 1 ,m2, . . . , m p }. The 
boundary which is the layer n = defines the preferential direction (the axis n). 

The Schrodinger equation can be rewritten as a recursion equation (in terms of the 
discrete-time n) 

1pn+l,m = {E — e n7m )lpn,m ~ V'n-l.m ~ ^ ^n,m'- (8) 

m' 

Summation over m' runs over the nearest neighbours of the site m. The on-site potentials 
e n m are independently and identically distributed. We assume hereafter existence of the 
two first moments, (e n , m ) = and (£^ m ) = cr 2 , where the parameter a characterizes the 
disorder level. 

Let us perform the Fourier transform: 

e n (k) =5>„, m e ikm , (9) 

m 

^n(k) = ^^, m e 4km . (10) 

m 

The relation for similar random quantities, , (e n) m) — 0, {s n ^ m e n ^ m i) = & 2 S ntn /5 mtm >, leads to 

(e n (k)) = 0, (11) 
( £n (k)<,(k')) = (27r)V 2 5 nX 5(k - k'). (12) 

As a result of the Fourier transform, the Schrodinger equation (jSJ) transforms into the equa- 
tion set for the mode dynamics, enumerated by the index k, 

VWi(k) = £(k)VUk) - Vw(k) - J S^(k - kOVn(ki). (13) 



Here 



p=D-l 



£(k) = E-2 c os(kj), 



(14) 



with fixed initial conditions ipo(k) and ipi(k). 

A comparison of eg. (1131) with its ID analog, eq.(j2J), demonstrates that the dynamics of the 
multi-dimensional system can be reduced to the dynamics of the multi-oscillatory system. 
The frequencies of these oscillators are defined by the relation uj(k) 2 = 2 — |£(k)|. That is, 
one can distinguish, as before, normal oscillators with |£(k)| < 2 and inverted oscillators 

n 

with |£(k)| > 2. As we have shown [16], the condition |£(k)| < 2 without perturbation 
corresponds to delocalized states inside the band; and contrary, |£(k)| > 2 corresponds to 
the solution outside the band. In other words, the connection between the classical oscillator 
type (normal or inverted) and the quantum-mechanical solution remains also for arbitrary 
dimensions. However, there is also an important difference between the ID and ND systems: 
in the former case the energy magnitude determines uniquely the oscillator type. The terms 
with a random force describe only the oscillator's stochastic self-interaction, due to its linear 
dependence, —e n q n in eq.fjSJ), the oscillator cannot be stopped, q = p = 0. 

In contrast, in the mult i- dimensional case the energy E no longer determines uniquely 
the system's state, the fixed initial conditions ipo(k) and ^i(k) define simultaneously the 
type and number of initially excited oscillators, among which can be found both normal 
and inverted oscillators. The integral in eq. (H3l) corresponding to a random force, describes 
now the stochastic interaction between oscillators. That is, knowledge of the ID system is 
not sufficient for the description of multi-dimensional systems; as we demonstrate below, 
fundamentally new effects arise here. 

Using the causality principle Q, [if]], the equation for the first momentum of a random 
amplitude (^(k)) is quite trivial: 



It is easy to see that for unstable modes, |£(k)| > 2, even the first moments are divergent, 
| (■?/>„ (k)) | — > oo, as n — > oo. Its analog in the classical phase space (p, q) corresponds to 
unbound trajectories. 

Since we associate the appearance of localized states with unbound trajectories for clas- 
sical oscillators, it is easy to formulate the necessary (but not sufficient) condition for the 



<^i(k)>=£(k)^ B (k)>-(^i(k)). 



(15) 



S 



appearance with disorder of the delocalized states: the initial conditions should correspond 
to excitation of normal oscillators only, i.e. the amplitudes ^o(k) and (k) are nonzero 



only for modes with |£(k)| < 2. This coincides with the statement [16| based on differ- 
ent ideas. Under this condition the dynamics of the first moments is bound for all modes, 
K'i/'nlk))! < oo as n — > oo. However, this condition is not sufficient, since the localized 
states as it was illustrated for the ID case, correspond in general to statistically unbound 
trajectories, (|^ n (k)| 2 ) — > oo as n — > oo. That is, the search for the sufficient condition for 
the excistence of the delocalized states is reduced to the solution of equations for the second 
— oftherando^plitudesflQ. 



C. Equations for second moments 

Divergence of the second moments is a typical diffusion behaviour; an observation of such 
a diffusion dynamics indicates directly the presence of the localized states, and vice versa. 
An easy criterium of a diffusion is the behaviour of the squared coordinate for all oscillators 

r dP~k 

Un = J ^(|^(k)| 2 ). (16) 

To detect the diffusion, it is sufficient to demonstrate divergence of the function U n as 
n — > oo. 

Let us define the second moments by the relations: 

x n (k) = (|^(k)| 2 ) = (^ n (k)<(k)), (17) 

y n (k) = i[(Vak)VCi(k) + C(k)VWi(k)>]. (18) 

One gets 

p^'„(k). (19) 



f.Qr 



Using the causality principle (see details in Ref. Il6|)) for eq. (jT3l) . one gets immediately 
for the non-zero average quantities: 

<|^ +1 (k)| 2 > = £ 2 (k)(|^(k)| 2 > + <|V„-i(k)| 2 > (20) 
-£(k)(Wk)Ci(k) + <(k)V>„_i(k)> 

f C APXr. WPb-o 

+ J J (27r)P (27r)P ^ k ~ k iX( k - k 2 ))(Vn(kx)V;(k 2 )). 
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Taking into account properties of random potentials, eq. (fT2l . one gets for n > 1: 

x n+1 (k) = £ 2 (k)x„(k) + x n _x(k) - 2£(k)y n (k) + a 2 f/ n . (21) 

Notice that the last term in eq. (l2~Tl) does not depend on k, i.e. the noise equally affects 
all modes. The noise intensity is described by cr 2 U n , where U n was defined in eq.(TT9l). 
Analogously, the complementary equation is derived 

y n+1 (k) + y n {k) = £(k)x„(k). (22) 

Let us perform now the Z-transform: 

*fek)=£^ (23) 

n=l 

nz , k) = f;^, (24) 

n=l Z 

oo „ 

n=l Z 

Takin into account eq. ffT9l . one gets easily 

f rl p k 

The Z-transform of eqs. (l2Tj) . (f22j) leads to the relations containing the initial conditions 

[z - z- 1 - C 2 (k)}X(z, k) + 2£(k)Y(z, k) = (27) 
xi(k) + 2T 1 x (k) + cx 2 £/(;z), 
(z + l)Y(z, k) - £(k)X(z, k) = Vl (k). (28) 



It is easy to find that 

(z-r 



[{z + lf/z - C 2 {k)}X{z, k) = \{z, k) + a 2 U(z), (29) 



{z + K 

where A(z,k) = 2£(k)yi(k) / (z + 1) + xo(k) + xi(k)/z. Taking into account the definition 
(|2T)j) . eq. (l2"91 is the integral equation, but with a simple structure. 

The initial conditions for moments are easily expressed through ?/>o(k), ^i(k). The trivial 
result is that if any mode with the wave vector ko enters the initial conditions (i.e. at least 
one of the two quantities ^o(k), 4>i(k) is non-zero), the quantity A (z, k ) is also non-zero for 
this mode. This will be used in the further analysis. 
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D. Localization operator 



Let us assume for the beginning that a = 0, i.e. there is no disorder in eq. (l29|) . and use 
the relevant solution X^°\z, k) for the calculation of the squared coordinate, eq. fl2"6T) . Thus, 
one gets 

U®(z) ~ fc±^ f ^ A( "' k) (30) 

For a 7^ the solution reads 

U(z) =H(z)U<®(z). (31) 

Here H(z) is the localization operator 

1 , 2 (z + l) f d p k 1 



H(z) " ° (z-l)J (27r)f [(z+l)Vz-£2(k)]' (32) 
Using the convolution property for the Z-transform Q, [3], one gets 

n 

Un = J2 U i° )h ^ ( 33 ) 
J=l 

where h n is the result of the inverse Z-transform of the localization operator H(z). 

As was mentioned, the necessary condition for the existence in the presence of disorder 
of the delocalized states is the presence at the initial time of normal (stable) modes only. 
In this case Un^ corresponds to the stable dynamics of the unperturbed problem and thus 
is bound in the index n: Un^ < oo as n — > oo. In its turn, the sufficient condition for the 
delocalized states is the absence of the (diffusion) U n divergence: U n < oo as n — > oo. As 
soon as U n — > oo, as n — > oo, this indicates the localized states. 

It is easy to notice that the convergence or divergence of U n is not dependent on the 
properties of the unperturbed solutions. The problem is reduced to a study of the asymptotic 
behaviour (n — ■> oo) of the h n coefficients in the linear transformation, eq. (l33l) 141 ]. This does 
not need calculation of the coefficients h n but, by the means of analytical methods, analysis 
of the localization operator H(z) as a function of the complex variable z. Therefore, the 
physical problem of the localized/delocalized states is reduced to the mathematical search 

l^l |"~ 



for the poles of the function H(z) of the complex variable 



, [15] . We do not go into details 



here; in particular, the diagrammic technique of the search for the poles was discussed in 



Appendix of Ref. 



id 
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In particular, it was shown Il5l ] that the localization operator H(z) is a non-analytic 
function of the complex variable z. The unit circle \z\ — 1 divides the complex plane into 
two analytic domains: the interior and exterior of the unit circle. The inverse Z-transform 
is quite generally defined via countour integrals in the complex plane 

This definition is only possible in an a naly tic domain and does not always represent a solution 
which can be physically interpreted [jjj, [jjj]. In this way, multiple solutions can result in 
the formal analysis of the problem. The first solution H + (z) describes the localized states. 
It is defined outside the unit circle and always exists. The second solution H_(z) describes 

re unit circle. The coexistence of the two solutions (if 



16| as the coexistence of two phases - sm insulating 



delocalized states. It is defined inside t 
any) was physically interpreted 
and a metallic one. 

Notice that the same result can be obtained using, instead of the total squared coordinate, 
eq. ffTB"]) . the total squared momentum as a criterium of the diffusion dynamics 



Vn = J ^(l^, + i(k)-^(k)| 2 ). (35) 

After simple transformations one gets 

V(z) = H(z)V {0) (z) } (36) 

whose structure is similar to that of eq. fl3~Tj) . the definition of the localization operator H(z) 
also retains, eq. (l32l . whereas V^°'(z) is the squared momentum of the unperturbed system. 
The expression for the localization operator H(z) has been derived by us earlier Q, 



151 ]. However, the nontraditional engineering language was used (such as input/output 
signals, filter function, etc). As we demonstrated above, the derivation can be considerably 
simplified, since in the system of interacting modes the asymptotic behaviour of the total 
squared coordinate U n serves as a natural indicator of the presence/absence of diffusion 
in the system's dynamics. Therefore, reformulation of the Schrodinger equation (quantum 
tight-binding model) in terms of the classical Hamiltonian map permits to retain the basic 

nrr 



definitions of the alternative approach 



16] , but opens the additional possibility of a 



new interpretation of the results obtained earlier, which we discuss in the next Section. 
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III. THE FPU PROBLEM VS THE AL PROBLEM 



A. Stability and thermalization 

It is generally believed that an increase in space dimension greatly increases the system's 
stability with respect to disorder. Thus, there is without doubt the presence of a metal- 
insulator transition in the 3D case. It is believed that the effect of statistical fluctuations 
changes the regime at D = 4 



29, 



30] ; no phase transitions are expected for D > 4. The 2D 



system marks the borderline between high and low dimensions 31 1 . However, all these con- 
elusions derive from the phenomenological scaling theory of localization [21[. The alternative 
point of view with different classification of high and low dimensions has been presented by 
us in Ref. 15|] . In general, this confirms the system stability in higher dimensions. 

An advantage of our interpretation of the AL in terms of classical dynamics of stochas- 
tically interacting oscillators is that it makes the statement on the relation between system 
stability and space dimension to be not so obvious. Indeed, for the energy range inside 
the old band \E\ < 2D, where normal modes certainly exist, simultaneously the inverted 
modes with |£(k)| > 2 are also always possible. Change of the space dimensions D affects 
the weights of these states, but their complete disappearance is impossible. That is, this 
problem can be resolved only by means of the exact analytical solution. 

On the other hand, the instability mechanism for multi-dimensional systems (D > 1) 
in terms of oscillators is quite obvious. Even assuming that at the beginning only normal 
modes are excited, stochastic interaction inavoidably excites also neighboring, in particular, 
inverted modes. In other words, thermalization of all modes takes place. It is important to 
stress that all modes - normal and inverted ones - contribute to the localization operator 
H(z), irrespective of which modes were excited in the beginning. The inverted mode dy- 
namics in the classical phase space (p, q) corresponds to unbound trajectories which do not 
correspond to the delocalized states. However, an immediate conclusion suggests itself that 
all solutions of the dynamical problem - independent of the space dimension and disorder 
level - correspond only to the localized states. Such a paradox conclusion demonstrates 
clearly that the problem under consideration is counter-intuitive. 

In our opinion, the AL paradox has much in common with the FPU problem [l|, [J . To 
show their close similarity, let us summarize here the main results. 
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Numerical simulations of a chain of harmonic oscillators coupled with a quadratic or cubic 
nonlinearity show that energy, initially placed in a low-frequency normal mode of the linear 
problem stay almost completely locked within a few neighbor modes (or quasi-modes Q]), 
instead of being distributed among all modes of the system. Recurrence of energy to the 
originally excited mode is also observed. The nonlinear effects are significant and cannot be 
neglected. 

Two alternative explanations of the FPU paradox were suggested {3]): the integrability 
of nonlinear equations and dynamical (deterministic) chaos. The second approach points to 
the existence of a stochasticity threshold in the FPU problem. If the nonlinearity is below 
a stochasticity threshold, the dynamics of the system remains similar to the one of the 
unperturbed system for large time scales. For a strong nonlinearity the overlap of nonlinear 
resonances leads to a strong dynamical chaos, destroying the FPU effect. Namely in this 
case the intuitive thermalization occurs. 

Therefore, in the standard FPU statement nonlinear effects play a key role being respon- 
sible for stochasticity: the motion of a nonlinear dynamical system even with few degrees 
of freedom can exhibit chaotic behavior. On the other hand, nonlinearity of the dynamics 
equations prevents their analytical analysis. A number of questions still remain open. In 
particular, the main results are obtained for systems with few degrees of freedom, the be- 
havior of the system is not known in the thermodynamic limit, when the number of degrees 
of freedom goes to infinity. 

In our interpretation of the Schrodinger equation in terms of the classical Hamiltonian 
map, stochasticity is introduced through on-site potentials, which are random variables. 
Therefore, there is no longer the need to solve nonlinear equations. For the Cauchy problem 
with fixed initial conditions, where the causality principle can be applied Q, [if]], this 
problem can be solved exactly analytically (see Section HlDj) . both for finite-size systems 
and its thermodynamic limit. Unlike the FPU problem with counter-intuitive solution, 
the classical interpretation of the AL problem is more complicated due to the unobvious 
dynamics of the inverted modes. 

We will consider below the delocalized states in the quantum mechanical problem as sta- 
tistically stable quasi-periodic motion, in terms of the classical Hamiltonian map, with ex- 
citations spreading only over a few neighbor modes (of the initially excited normal modes). 
Such a dynamics is bound. The localized states are interpreted respectively as diffusion 
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dynamics with statistically unbound trajectories. In the FPU problem this means thermal- 
ization with an excitation spreading over all modes. 

B. Stochasticity threshold vs disorder threshold 

As was noticed in Ref. [si], the existence of a stochasticity threshold in the FPU is the 
nonlinear effect. For a strong nonlinearity the overlap of nonlinear resonances leads to a 
dynamical chaos, destroying the FPU recurrence. As a result, fast convergence to thermal 
equilibrium arises. The problem of the nonlinearity prevents its detailed analytical study. 
The stochasticity threshold depends on the type of nonlinearity, space dimension; the ther- 
modynamic limit is also unclear. Compared to this situation, the dynamical interpretation 
of the Schrodinger equation has obvious advantages permitting an exact analytical solu- 
tion. It is expected that in the stochastic AL problem the relevant stochasticity threshold 
transforms into the disorder threshold. 

When comparing the two problems, one has to keep in mind that the AL was formulated 
as a statistical problem with an ensemble of random potential realizations. Thus, it is 
convenient to reformulate the FPU as a statistical problem. Introducing the stochasticity 
parameter K (nonlinearity coefficient [26]), the type of dynamical trajectories is defined by 
the initial conditions. In other words, in the nonlinear dynamic problem the stochasticity 
parameter is not a unique factor determining the dynamics. For example, the behavior of 
the FPU system depends strongly on whether low- or high-frequency modes are initially 
excited The number of excited modes seems also to be important for the dynamics. 

As is well known, in such problems the phase space is divided into regions with qual- 
itatively different types of motion, and these regions are separated by barriers. If initial 
conditions were chosen in the region corresponding to a stable quasi-periodic motion, this 
dynamics corresponds to the recurrent behavior as in the FPU experiment and is classified 
as the dynamics below the stochasticity threshold. In another region an instability of mo- 
tion exists for a wide range of the initial conditions (the dynamics above the stochasticity 
threshold). The value of the stochasticity parameter K determines the borders of these 
regions. 

Let us now define some domain in phase space and an ensemble of the initial conditions 
therein. If for a given K a whole domain chosen falls into some region, all trajectories in 
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this ensemble reveal the same dynamics - stable quasi-periodic motion or dynamical chaos 
(single phase domain) - otherwise an ensemble reveals trajectories of different kinds (two- 
phase domain). In the latter case trivial co-existence of the two phases (or two dynamics, 
in a statistical sense) takes place. Change of the stochasticity parameter K in a given 
ensemble leads to the phase transitions. It is easy to show that in our statistical problem 
phase transitions from single-phase to hetero-phase (phase co-existence) dominate, and vice 
versa, since these transitions correspond to the boundary motion of the regions. Transitions 
from one single-phase system to another single-phase are also possible since variation in the 
K parameter in the phase space can induce spontaneous creation of new regions with a 
different dynamics. 

The Schrodinger equation in terms of the classical Hamiltonian map does not depend 
so strongly on the initial conditions. In fact, these determine only the quantity defined 
by eq. (l50l) . whereas existence of localized/delocalized states is defined entirely by the lo- 
calization operator H(z), eq.( l32l) . which does not depend at all on the initial conditions 
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161 ] . The physical reason for this is obvious. The ensemble of different trajectories 



in the AL can be created already for fixed initial conditions by means of the ensemble of 
random potential realizations. For these realizations, even qualitatively different trajectories 
in the classical phase space occur. There is no reason to believe that the disorder parameter 
a defines uniquely the trajectory types. It can be assumed that some random potentials 
(called coherent realization) correspond to the solutions ip n (\i) close to average values (0 n (k)) 
(the delocalization regime). For these realizations excitation of neighbor modes also occurs, 
however, at the next discrete-time steps n these modes return to an equilibrium position, 
ip n (k) = (the recurrent behavior). In other words, only their virtual dynamics around 
equilibrium with limited amplitude takes place. For other potential realizations (the lo- 
calization regime) thermalization occurs and the recurrent behavior disappears. The main 
physical question - whether these two regimes have comparable statistical contributions - 
needs drawing of the phase diagram 

Since the detailed AL study in the thermodynamic limit for an arbitrary space dimension 



D has already been published by us ,15], we restrict ourselves here to the interpretation of the 
result from the point of view of the disorder threshold. As was shown [15fl . for D > D = 4 
the spectrum of wave function fluctuations changes, as well as the convergence of the integral 
in eq. (|32l) . The spatial dimension Dq = 4 was also discussed in Ref.Q, Q , however, the 
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conclusion was drawn therein that this is an upper critical dimension for localization, so 
no phase transitions are expected for D > 4. In other words, no localized states can exist 
here. Since appearance of the localized states in the dynamical version of the Schrodinger 
equation ([TBI means thermalization (i.e. energy transfer from initially excited modes to all 
other modes), the analog of the above-mentioned result [30] in the FPU is a very strong 
statement on the existence of entirely recurrent states and no convergence to the thermal 
equilibrium in dynamical systems for D > 4. That is, the Gibbs statistics would not be 
applicable for high-dimensional systems. However, this statement is not based on an exact 
solutions and looks very suspicious to us. Thus, let us consider the alternative idea. 

According to our Ref. 15], for D > D = 4 the problem is fundamentally simplified: 
for all energies inside the old band \E\ < 2D (where only the delocalized states exist) 
there is a disorder threshold ctq(E). There is no thermalization for a < <jq(E), but only 
delocalized states, characterized by the (formal) Lyapunov exponent 7 = 0. For a > cr (E) 
thermalization occurs, all states are exponentially localized with a certain j(a,E). As 
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16] , the Lyapunov exponent 7 can be treated in the AL problem as the long- 



range order parameter: the two different phases reveal different 7 values, 7 = (delocalized 
states) and 7 7^ (localized states, respectively). A sharp separation in the phase diagram of 
the localized and delocalized states recalls the second-order phase transition, where phases 
cannot coexist. However, the phase transition here is not of the second-order, since 7 reveals 
step-like changes at the disorder threshold. As the space dimension D increases, the system 
stability increases with respect to disorder, but thermalization in principle cannot completely 
disappears. 

For the case 2 < D < 3 (low-dimensionial case 15]) the problem becomes more compli- 
cated; the energy threshold Eq = Eq(D) also arises. The delocalized states disappear even 
under infinitesimal disorder on the boundaries of the old band, E < \E\ < 2D. In other 
words, only thermalization is possible in this energy range. There is no analog of the FPU 
recurrence here. In the region \E\ < Eq the disorder threshold <Jq{E) again appears; e.g. 
for D = 2 E = 2 and a (E) = 2(1 - a 2 /^ 2 ) 1 / 4 [14]), but its sense changes. If disorder a 
exceeds the threshold, only the localized states exist (full thermalization). As a < a (E), 
the disorder parameter a no longer separates uniquelly the localized and delocalized states, 

. In other words, in the ensemble of differ- 



since phase co-existence takes place 
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ent disorder potential realizations in this region both localized and delocalized solution can 
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arise with comparable probabilities. As a result of phase co-existence, phase transition in 
low-dimensional systems should be considered as a first-order transition. As was said above, 
this non-trivial result has a direct analogy in the statistical version of the FPU problem. 

In the ID case the analogy between the AL and FPU problems disappears; the 
Schrodinger equation corresponds here to a single oscillator only and thermalization has 
no longer a meaning. Notice, however, that historically the FPU problem arose during the 
numerical simulations of a one-dimensional chain of oscillators Q, S|. The localization in ID 
means diffusion excitation in one degree of freedom with specific disorder effects. This is 
not surprising: when presenting the Schrodinger equation in the form of eq.( fT3|) . one of the 
spatial coordinates in the D-dimensional system was interpreted as a discrete time variable, 
respectively, the Hamiltonian dynamics is treated for oscillator systems in a p = D — 1 space 
dimension. That is, the AL problem in a space dimension D should be compared with the 
FPU with a lower space dimension, p = D — 1. 

C. FPU recurrence vs delocalization 

For the sake of illustration, we restrict ourself to the case of the band centre E = for 
2D system. For this particular case for the delocalized states (H(z) = H-(z)) we calculated 
analytically [l4j the inverse Z-transform, eq. (IMl) . and found the coefficients h n in eq. (l33l . 
Eq. (l33l) can be transformed, respectively: 



Eq. (l3"TI) is valid only for a < 2, provided 2sin(0) = a [141 ] . 

It is easy to notice that the solution for a disordered system derives from that for the 
ordered system with the help of the sinusoidal modulation. Taking into account that the 
full squared coordinate Un is a sum of oscillating quantities (normal modes), solution for 
the perturbed problem U n is a quasi-periodic motion. That is, one observes a direct analog 
of the recurrent behaviour in the FPU systems. 

As was mentioned above, the 2D case for the AL corresponds to the ID dynamics. On 
the other hand, the FPU paradox was established for the first time namely for the one- 
dimensional system , and the result was considered as counter-intuitive. This is why 
our exact analytical study of the AL in the 2D case and the statement on the existence of the 




(37) 
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delocalized states [14| proves once more a close similarity between the AL and FPU. Indeed, 
if the recurrent behavior is possible for the ID FPU and this result is generally accepted as 
reliable, there are no strong grounds to reject our conclusion on the delocalized states in the 
2D AL, since these states are analogous to the ID quasi-periodic motion in FPU. 

IV. CONCLUSION 

The Anderson localization and Fermi-Pasta-Ulam problems are very complementary; each 
non-trivial and, as a rule, counter-intuitive result in one problem has its analog in another 
problem. In particular, a stable quasi-periodic motion (the recurrent behavior) in the FPU 
corresponds to the existence of the delocalized states in the AL problem. In contrast, the 
thermalization effect in FPU has an analog in the localized states in the AL. In general, we 
have shown here that a deep analogy between these two problems is possible, when we treat 
the quantum tight-binding model in terms of a classical Hamiltonian map. 

This work was partially supported by the Latvian Council of Science. Author is greatly 
indebted to W. von Niessen for detailed discussions of the paper. 
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